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According to a recent proposal, the so-called Barbero-Immirzi parameter of Loop Quantum Grav- 
ity can be fixed, using Bohr's correspondence principle, from a knowledge of highly-damped black 
hole oscillation frequencies. Such frequencies are rather difficult to compute, even for Schwarzschild 
black holes. However, it is now quite likely that they may provide a fundamental link between 
classical general relativity and quantum theories of gravity. Here we carry out the first numerical 
computation of very highly damped quasinormal modes (QNM's) for charged and rotating black 
holes. In the Reissner-Nordstrom case QNM frequencies and damping times show an oscillatory 
behaviour as a function of charge. The oscillations become faster as the mode order increases. At 
fixed mode order, QNM's describe spirals in the complex plane as the charge is increased, tending 

£f*) • towards a well defined limit as the hole becomes extremal. Kerr QNM's have a similar oscillatory 

behaviour when the angular index m — 0. For I = m = 2 the real part of Kerr QNM frequencies 
tends to 2f2, f2 being the angular velocity of the black hole horizon, while the asymptotic spacing 

£SJ ■ of the imaginary parts is given by 2-kTh- 

g \ 

' 1. INTRODUCTION 

m ■ 

Quasinormal modes (QNM's) play a fundamental role in black hole physics. They are known to "carry the fin- 
gerprints" of a black hole, since their frequencies only depend on fundamental black hole parameters such as mass, 
charge and angular momentum. They determine the late-time evolution of fields in the black hole exterior. Even 
more importantly, they may play a fundamental role in the newborn field of gravitational wave astrophysics; indeed, 
numerical simulations of stellar collapse and black hole collisions have shown that in the final stage of such processes 
("ringdown") QNM's dominate the black hole response to any kind of perturbation. Therefore QNM's have been 
Qf~^ . extensively studied for more than thirty years (for comprehensive reviews see 0,0) with the aim to shed light on the 
f^*) ' strong-field behaviour of classical general relativity. 

Research was mainly focused on modes having small imaginary part for two main reasons: first of all, weakly 
damped modes are expected to be dominant in gravitational wave radiation; in the second place, numerical methods 
used to compute QNM's generally run into trouble when the modes' imaginary part grows, i.e., when the damping 
is high. Early investigations of highly-damped modes had motivations which were radically different from those of 
O .' the present paper. Highly-damped modes were seen as a good benchmark for the reliability of numerical methods. 
Their study could provide hints (if not a formal proof) to whether QNM's are or not infinite in number. Finally, 
1 there was hope that their study could lead to a better understanding of the long-standing issue of mode completeness. 
J> \ Applications of different methods yielded at first puzzling and contradictory results: WKB methods predicted that the 
'kJi ' asymptotic real part of the frequency, lur, should vanish for highly damped modes [3(, while the continued- fraction 
methods developed by Leaver |4| seemed to suggest that ojr should be finite. An improvement of the continued 
fraction technique devised by Nollert 5] finally showed that ujr is indeed finite and determined its value; these results 
were confirmed by Andersson 6] building on an improved WKB- type technique, the phase integral method, previously 
developed in 0. 

More recently QNM's, which are essentially related to the classical dynamical properties of a black hole, have 
become a subject of great interest for the quantum gravity community. This interest stems essentially from Hod's 
proposal to apply Bohr's correspondence principle to black hole physics. By the Bekenstein-Hawking formula the 
surface area of a black hole is nothing but its entropy. In a quantum theory of gravity the surface area should have a 
discrete spectrum, and the eigenvalues of this spectrum are likely to be uniformly spaced. In order to give a prediction 
on the area spacing, Hod observed that the real parts of the asymptotic (highly damped) quasinormal frequencies of 
a Schwarzschild black hole of mass M, as numerically computed by Nollert |5(, can be written as 



u>R = T H ln3. (1) 

where Tjj = (87rM) _1 is the black hole Hawking temperature (here and in the following we use units such that 
c = G = 1). He then exploited Bohr's correspondence principle, requiring that transition frequencies at large quantum 
numbers should equal classical oscillation frequencies, to predict the spacing in the area spectrum for a Schwarzschild 
black hole. It is worth stressing that in this quantum gravity context, as opposed to the study of QNM's in the 
context of gravitational wave emission, relevant modes are those for which the imaginary part tends to infinity - that 
is, modes which damp infinitely fast and do not radiate at all. 
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We point out that a dynamical interpretation of the Hawking effect through a semiclassical treatment of the 
quantum- mechanical uncertainty associated to QNM oscillations was proposed long ago by York . York's dynamical 
treatment could successfully be used to calculate the temperature and entropy of the hole. However, his proposal to 
relate classical black hole oscillations to the hole's quantum-mechanical behaviour is essentially different from Hod's, 
in that the dynamical effects considered in are dominated by slowly damped QNM's. 

Following Hod's suggestion, Bohr's correspondence principle has recently been used by Dreyer to fix a free parameter 
(the so-called Barbero-Immirzi parameter) appearing in Loop Quantum Gravity |10| . Supposing that transitions of a 
quantum black hole are characterized by the appearance or disappearance of a puncture with lowest possible spin j m i n , 
and that changes AM in the black hole mass M, corresponding to such a transition, are related to the asymptotic 
frequency by 

AM = fiio R , (2) 

Dreyer found that Loop Quantum Gravity gives a correct prediction for the Bekenstein-Hawking entropy if j m in = 1 , 
consequently fixing the Barbero-Immirzi parameter. Motivated by the occurrence of an integer value for j m in, Dreyer 
went on to suggest that the gau ge g roup of Loop Quantum Gravity should be SO(3), and not SU(2). Such a proposal 
has recently been questioned in and |T^ . In the latter paper, Hod's proposal has also been used as an argument 
in favour of an equidistant black hole area spectrum. Evidence that the earlier formula for black hole entropy in Loop 
Quantum Gravity still holds when j min = 1 has been presented in p| . 

Numerical computations of gravitational QNM's for higher-dimensional black holes are still lacking. However, 
Kunstatter [l4j recently generalized Dreyer's argument to give a prediction for the area spacing (and for the asymptotic 
oscillation frequency) of d-dimensional black holes. 

When Hod made his original proposal, the fact that lor is proportional to In 3 was just a curious numerical 
coincidence. Further support to the aforementioned arguments came from the analytical proof by Motl |15| that 
highly damped QNM frequencies are indeed proportional to In 3. More precisely, choosing units such that 2M = 1, 
high-n mode frequencies satisfy the relation 

uSchw ^ ln3 i(n-V2) 1/2 

47T 2 

Such an expression has recently been confirmed through a different analytical approach in |16j . where it has also 
been generalized to higher-dimensional black holes and to four-dimensional Reissner-Nordstrom (RN) black holes. The 
WKB approach used in ^(| has later been exploited to analytically compute reflexion and transmission coefficients of 
multidimensional Schwarzschild black holes (and of four-dimensional RN black holes) in the limit of large imaginary 
frequencies [l7j . 

There are many important reasons to try to understand the behaviour of highly-damped QNM's for general black 
holes. Let us consider charged, rotating black holes, having angular momentum per unit mass a = J/M and charge 
Q. The black hole's (event and inner) horizons are given in terms of the black hole parameters by r± = M ± 
\J M 2 — a 2 — Q 2 . The hole's temperature Tjj = {r+ — r-)/A, where A — 47r(2Afr + — Q 2 ) is the hole's surface area, 
related to its entropy S by the relation A = 5/4. Let us introduce the so-called "angular velocity of the horizon" 
n = Aira/A, and $ = AnQr + /A. Applying the first law of black hole thermodynamics, 

AM = T H AS + QAJ + $AQ, (4) 

but dropping without justification the AQ term, and assuming that the formula for the area spectrum derived for a 
Schwarzschild black hole still holds in this case, Hod conjectured @ that the real parts of the asymptotic frequencies 
for charged and rotating black holes are given by: 

lor = ujji = Th In 3 + mO, (5) 

where m is the azimuthal eigenvalue of the field. In particular, such a conjecture implies that QNM's of extremal RN 
black holes would have a vanishing asymptotic real part. Unfortunately, numerical studies of asymptotic frequencies of 
charged and rotating black holes have been lacking until now. Hod |lg recently used the most systematic exploration 
of Kerr black hole QNM's, which was carried out a few years ago by Onozawa [l9|. to lend qualitative support to 
formula (J5J. However, in the following we will extend Onozawa's numerical calculations to larger imaginary part, 
showing that the use of low-order frequencies to deduce the asymptotic behaviour as u>i — > oo is rather questionable. 

Motl and Neitzke |16| recently obtained an analytic formula for the asymptotic frequencies of scalar and 
electromagnetic-gravitational perturbations of a RN black hole: 

e PuJ + 2 + 3e _/3/w = 0. (6) 
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For computational convenience, the authors fixed their units in a somewhat unconventional way: they introduced a 
parameter k related to the black hole charge and mass by Q/M = 2Vfc/(l + k), so that (3 = An /(I — k) = 1/Th is 
the inverse black hole Hawking temperature and (3i = —k 2 (3 is the inverse Hawking temperature of the inner horizon. 
However, some features of their result are particularly puzzling: 

1) As emphasized (and partially justified with plausibility arguments) in [Tiil. Il7|. the predicted asymptotic RN 
quasinormal frequencies do not reduce to the Schwarzschild limit as the black hole charge Q tends to zero; 

2) Quasinormal frequencies of a charged black hole, according to formula ©, depend not only on the black hole's 
Hawking temperature, but also on the Hawking temperature of the (causally disconnected) inner horizon; 

3) The authors suggest that, should the black hole mass and charge acquire a small imaginary part (which in their 
words "may not be an unreasonable thing to do", since "the black hole eventually evaporates"), their asymptotic 
RN frequencies would be proportional to In 2. This is in stark contrast with the Schwarzschild results: should 
this be true, Dreyer's argument could be used to infer that the gauge group of Loop Quantum Gravity is SU(2); 

4) The result does not seem to agree with the conjectured behaviour predicted by formula if5|. 

Therefore, their result cannot be considered conclusive, and there are many issues to clarify. Furthermore, if 
asymptotic frequencies for "generic" black holes depend on the hole's charge and angular momentum, a relevant 
question is: how should arguments based on Bohr's correspondence principle be modified? A similar question was 
recently raised by Cardoso and Lemos |20j | . They studied the asymptotic spectrum of Schwarzschild black holes in de 
Sitter spacetimes and found that, when the black hole radius is comparable to the cosmological radius, the asymptotic 
spectrum depends not only on the hole's parameters, but also on the angular separation index I. 

For all these reasons, a numerical computation of highly-damped QNM's is now needed more than ever. Such a 
numerical computation is technically challenging: even though the first overtones of the QNM spectrum have now 
been studied for more than 15 years, we present here the first computation of this kind for RN and Kerr black holes. 

We will present numerical results that generally support the calculations carried out in 0, H3, . Indeed, a 
testable prediction of those analytical derivations is that asymptotic frequencies for scalar perturbations should have 
the same value as gravitational frequencies; however, to our knowledge, asymptotic scalar modes have never been 
shown in the published literature. We will extend Nollert's calculation to scalar modes, confirming the analytical 
prediction. Furthermore, our scalar mode calculation gives useful hints on leading order corrections to the asymptotic 
frequencies of Schwarzschild black holes. 

The plan of the paper is as follows. In section [21 we describe our numerical method, extending Nollert's technique 
to RN and Kerr black holes. In section [3] we show our numerical results, and put forward some conjectures on their 
implications for the asymptotic behaviour of the modes. The conclusions and a discussion follow. 



2. NUMERICAL METHOD 

A first comprehensive analysis of the QNM spectrum of RN black holes was first carried out by Gunter [2l| : 
then Kokkotas and Schutz [23 verified and extended his results using numerical integrations and WKB methods. 
Unfortunately, the standard WKB techniques cannot be applied to our case, since they become inaccurate in estimating 
the modes' real part as the imaginary part increases, unless one resorts to more sophisticated phase-integral methods 
• The first few modes of the Schwarzschild and Kerr black holes were studied by Leaver using a continued fraction 
technique which was then extended to the RN case |2^]. This technique is generally rather accurate for modes 
having u>i ~ u>r, but it eventually loses accuracy when u>i 3> u>r. The error is essentially introduced by a truncation 
of the power-series solution to the radial equation at some large N. For large N, Leaver has shown that this error 
can be written as an integral, which is rapidly convergent near the lower quasinormal frequencies (where \u>r\ > \u>i\ 
and \uj\ ~ 1); however, the convergence becomes slower as the overtone index increases [23. That's why the problem 
of computing high-order overtones is so numerically challenging. An improvement of the continued fraction method 
was used by Nollert 5] to find the asymptotic behaviour of the QNM frequencies for Schwarzschild black holes. In 
the following sections we show how to generalize Nollert's method to the charged and rotating cases. 
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2.1. Reissner-Nordstrom black holes 



In this section and the following we briefly describe the computational procedure we used. More details are given, 
eg., in 0, H, E3- Let us introduce a tortoise coordinate r*, defined in the usual way by the relation 

£ = i P) 

ar* r A 

where A = r 2 — 2Mr + Q 2 ; after a separation of the angular dependence and a Fourier decomposition, axial electro- 
magnetic and gravitational perturbations of a RN metric are described by a couple of wave equations: 

^ + ^)z r = Vr Zr . (8) 

Polar perturbations can be obtained from the axial ones through a Chandrasekhar transformation j24|. In the limit 
Q = 0, the potentials V{~ and V 2 describe, respectively, purely electromagnetic and axial gravitational perturbations 
of a Schwarzschild black hole. Now, the radial equations for the perturbations can be solved using a series expansion 
around some suitably chosen point. The coefficients a n of the expansion are then determined by a recursion relation. 
For Schwarzschild black holes the recursion relation has three terms, i.e., it is of the form: 

a ai + f3 a = 0, (9) 
aia 2 + p\a\ + 7ia = 0, 

a n a n+ i + f3 n a n + 7„a„_i + <5„a„_ 2 = 0, n = 2, 3, . . . 

where the recursion coefficients a n , f3 n and j n are functions of the frequency u> and of I (if we fix units such that 
2M=1). 

For RN black holes we actually have a four-term recursion relation (whose coefficients also depend on the charge 
Q), but we can reduce it to the previous form using a Gaussian elimination step. It turns out that the QNM boundary 
conditions are satisfied when the following continued-fraction condition on the recursion coefficients holds: 

= /? -^f^... (10) 

Pl— P2 — 

The n-th quasinormal frequency is (numerically) the most stable root of the n-th inversion of the continued-fraction 
relation QlOfl. i.e., it is the root of 

Pn r, o ■ ■ ■ n ~ y^^> 

Pn-\— Pn-2~ PO 

= - R a •■■ (n = l,2,...). 

Pn+l— Pri+2 — 

The infinite continued fraction appearing in equation l|llf) can be summed "bottom to top" starting from some large 
truncation index N. Nollert has shown that the convergence of the procedure improves if such a sum is started using 
a wise choice for the value of the "rest" of the continued fraction, Rn- This rest can be obtained solving the equation 

Rn = 1N+1 R , (12) 

and assuming that 

oo 

Rn = J2 CkN ~ k/2 - ( 13 ) 

fe=0 

If we introduce p = — iw, and denote by r+ the coordinate radius of the black hole's event horizon, it turns out that 
the first few coefficients in the series are Co = — 1, C\ = ±^/2p(2r + — 1), Ci — (3/4 — 2pr + ). 

As first shown in using WKB techniques and then confirmed using different numerical methods in j2|| , even for 
modes with moderate imaginary parts, RN QNM frequencies show a very peculiar behaviour as the charge increases. 
Our numerical results agree with those shown in [26|. Furthermore, we have checked that our method yields Nollert 's 
asymptotic frequencies in the Schwarzschild limit. Due to convergence reasons, calculations become more and more 
computationally intensive as the charge is increased. Indeed, because of the merging of the inner and outer horizons, 
maximally charged black holes yield radial equations for the perturbation variables which have a different singularity 
structure, and deserve a special treatment [23. We plan to investigate the behaviour of highly damped QNM's for 
maximally charged black holes in the future. 
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2.2. Kerr black holes 



As we did for the charged black hole case, here we only briefly sketch our computational procedure, referring the 
reader to 0, 0, 0] for more details. In the Kerr case, the perturbation problem reduces to a couple of differential 
equations - one for the angular part of the perturbations, and the other for the radial part. In Boyer-Lindquist 
coordinates, defining u = cos#, the angular equation reads 



[(1 



(auiu) 2 — 2ausu + s + A[ r - 



(m + sit) 2 



Sir, 



0. 



(14) 



and the radial one is 



Ai? Zm , rr + {s + l)(2r - l)Rim, r + V(r)R lm = 0, (15) 

where 

V(r) = {[(r 2 + a 2 fuj 2 - 2amusr + a 2 m 2 + is(am(2r - 1) - uj(r 2 - a 2 ))] A -1 

+ [2isur - a 2 u 2 - A lm ] } . (16) 

In writing down these equation we have adopted, consistently with what we did for the RN case, Leaver's con- 
ventions. In particular, we chose units such that 2M = 1. The parameter s = 0, — 1, —2 for scalar, electromagnetic 
and gravitational perturbations respectively, a is the Kerr rotation parameter (0 < a < 1/2), and Ai m is an angular 
separation constant. In the Schwarzschild limit the angular separation constant can be determined analytically, and 
is given by the relation Ai m = 1(1 + 1) — s(s + 1). 

Boundary conditions for each equation translate into a couple of three-term continued fraction relations of the 
form IjlOfl . Finding QNM frequencies is now a two-step procedure: for assigned values of a, £, m and uj, first find the 
angular separation constant Ai m (uj) looking for zeros of the angular continued fraction; then replace the corresponding 
eigenvalue into the radial continued fraction, and look for its zeros as a function of w. In principle, the convergence 
of the procedure for modes with large imaginary parts can be improved, as described earlier, by a wise choice of the 
rest, Rn, of the radial continued fraction. Expanding this rest as in formula (|13|) and introducing b = yl — 4a 2 , we 
get for the first few coefficients: Co = — 1, Ci = ±*/2pB, Gi = [3/4 — p(b + 1) — s]. 

As for the RN case, we have checked that our results agree with those shown in 0,^1 f° r small values of u>i, and 
that we get Nollert's asymptotic quasinormal frequencies in the non rotating limit. 



3. RESULTS 



As a first step in the computation of asymptotic modes for charged and rotating black holes we have checked that 
our numerical methods reproduce known results in the Schwarzschild limit |28j |. In particular, we have verified the 
asymptotic behaviour found by Nollert [j| . His main result was that the real part of the asymptotic QNM frequencies 
corresponding to gravitational perturbations can be well fitted by a relation of the form 

u>fl = w 0O + ^. (17) 

The leading-order fitting coefficient is independent of I and given by = 0.0874247, consistently (within numerical 
accuracy) with the analytical formula J2J). Corrections of order ~ n -1 / 2 , however, are i-dependent. Furthermore, we 
will see in a moment that they also depend on the spin s of the perturbing field, and that's why we denoted them 
by A s .;. For gravitational perturbations (s = —2) Nollert found (and we verified to the same level of accuracy) that 
A_ 2 ,2 = 0.4850, A_ 2 ,3 = 1-067, A_ 2 ,6 = 3.97. 

Now, an important testable prediction of the recent derivations of the "In 3" asymptotic behaviour is that scalar 
black hole perturbations should lead to the same asymptotic QNM frequency, at leading order in an expansion in 
powers of n -1 / 2 pH llii. Il7| . To check this prediction we have extended our numerical calculations to s = 0. Using 
the procedure described in we found that our scalar QNM data are well fitted by formula l|17[) . As n — > 00, the 
asymptotic frequency for scalar modes is again given by Woo = 0.0874247, consistently with the analytic calculation. 
What changes are the numerical values of the leading-order correction coefficients Ao,;. Namely, we find: Aq.o = 0.0970, 
A ,i = 0.679, Ao, 2 = 1-85. 
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Neitzke recently suggested that leading order corrections to the asymptotic frequency should be proportional 
to [(s 2 — 1) — 31(1 + 1)]. As an interesting by-product of our calculation, we found that our numerical values for the 
\ s /s are consistent with Neitzke's conjecture; that the proportionality constant is dependent on s, 

\ s ,i = k s [(s 2 - 1) - 31(1 + I)}; (18) 

and finally, we determined the proportionality constants k s to be given by fc_2 = —0.0323, ko = —0.0970. Recently 
Maassen van den Brink |29j derived fc_2 analytically, finding 

k - 2 = * -°- 0323356 - ( 19 ) 

Notice that, within our numerical accuracy, ko = 3fc_2- This prediction may be helpful to analytically determine ko. 
We shall say more about the importance of leading-order corrections to the frequencies in the following; now we turn 
to a discussion of our numerical results for RN and Kerr black holes. 



3.1. Reissner-Nordstrom black holes 



The numerical behaviour of the first few overtones of a RN black hole was studied numerically by Andersson and 
Onozawa 26], who unveiled a very peculiar behaviour as the modes' imaginary part increases. Our numerical codes are 
in excellent agreement with their results. In comparing with their paper note, however, that Andersson and Onozawa 
count modes starting from n — 0, while we label the fundamental mode by n — 1, following Leaver (see table 1 in 
0). The results we display refer to perturbations reducing to pure gravitational perturbations of Schwarzschild in the 
uncharged limit. The trajectories described by the modes in the complex-cj plane first show "closed loops", as in the 
top left panel of figure^ Then they get a spiral-like shape, moving out of their Schwarzschild value and "looping in" 
towards some limiting frequency as Q tends to the extremal value. This kind of behaviour is shown in the top right 
panel of figure^ We have observed that such a spiralling behaviour sets in for larger values of the modes' imaginary 
part (i.e., larger values of n) as the angular index I increases. In other words, increasing I for a given value of the 
mode index n has the effect of "unwinding" the spirals, as can be observed in the two panels in the second row of 
figure^ However, for each / the spiralling behaviour is eventually observed when n is large enough; a typical example 
for I = 6 is shown in the last plot of figure ^ 

Maybe a clearer picture of the modes' behaviour can be obtained looking separately at the real and imaginary 
parts of the mode frequencies as function of charge. Let us focus first on the real parts. The corresponding numerical 
results are shown in figure |3 together with the predictions of the analytic formula JSJ. It's quite apparent that, as 
the mode order grows, the oscillating behaviour as a function of charge start earlier and earlier. As n increases, the 
oscillations become faster, the convergence of the continued fraction method slower, and the required computing time 
gets longer. Therefore, when the imaginary part increases it becomes more difficult to follow the roots numerically 
as we approach the extremal value Q = 1/2. That's why our data for large values of n do not cover the whole range 
of allowed values for Q. Despite these difficulties, we have many reasons to trust our numerics. We have carefully 
checked our results, using first double and then quadrupole precision in our Fortran codes (indeed, as n increases, 
we can obtain results for large values of the charge only using quadrupole precision). As we have shown earlier, our 
frequencies accurately reproduce Nollert's results in the Schwarzschild limit, so our numerics can be trusted for small 
values of Q. Furthermore, the predictions of the analytic formula exactly overlap with the oscillations we observe 
for large values of the charge. Not only this does give support to the asymptotic formula 0]; ^ also gives us faith 
that the numerics are meaningful also for large charge, where we don't have any published results to confirm our 
predictions. 

Similar considerations apply to the imaginary parts. We present a few plots illustrating the general trend for the 
imaginary parts in figure [3] showing again excellent agreement with the asymptotic formula © as n increases. Once 
again, the analytic formula shows deviations from our "exact" numerical results only for small values of the charge, 
probably indicating that corrections of order n" 1 / 2 should be taken into account for small values of Q. When we 
look at the mode trajectories in the complex-w plane, the increasingly oscillating behaviour of the real and imaginary 
parts means that the number of "spirals" described by the mode before reaching the extremal value increases roughly 
as the mode order n. 

Can we deduce something from the agreement of our numerical results with formula 10 at large values of the 
charge? It would be extremely interesting to draw consequences on the extremal RN case, for various reasons. First 
of all, the QNM spectrum for extremal RN black holes is characterized by an isospectrality between electromagnetic 
and gravitational perturbations, which has been motivated in |30j as a manifestation of supersymmetry. Furthermore, 
topological arguments have been used to show that the entropy-area relation breaks down for extreme QNM's [3l| . 
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Therefore, we believe that some caution is required in claiming that the connection between QNM's and the area 
spectrum is still valid for extreme black holes, as recently advocated in |3^ |. These problems may be connected with 
our recent finding that extremally charged black holes in a (non-asymptotically flat) anti-de Sitter spacetime could 
be marginally unstable [3j|. 

Our numerics seem to indicate that we can trust formula Jfjjl in the large-charge regime. Then a very interesting 
conclusion follows 17]: the real part of the frequency for extremal RN black holes coincides with the Schwarzschild 
value, i.e. 

„ N ^ln3 ag q^ 1/2 ^ (20) 
An 

However, it is quite difficult to check this prediction numerically. In the extremal RN case, due to the coalescence of 
the inner and outer horizons, the singularity structure of the radial perturbation equations changes. Therefore one 
has to apply a different (and slightly more involved) procedure, which has been described in |27| . We have tried to 
apply that procedure to get highly damped QNM's. Apparently, the extremal RN QNM's show a behaviour which is 
rather similar to that of Schwarzschild QNM's (see, e.g., figure 1 in Q): they have finite real part for small values of 
luj, approach the pure- imaginary axis as a->j is increased, and then the real part increases again. This would support 
the predictions of the asymptotic formula (jfJJl . Unfortunately, we have not yet managed to get stable numerical results 
for large values of uij. We plan to improve our codes and obtain more numerically stable results in the future. 

If, supported by the agreement between our numerics and the analytical prediction, we assume that formula i|20|) 
holds, an interesting result emerges: Hod's conjecture J5| is incompatible with the "truncated" version of the first law 
of black hole mechanics l@J obtained dropping the AQ term Q- Indeed, formula (JHJ) predicts that the real part of the 
asymptotic frequencies in the extremal case should be zero. This does not imply that Hod's conjecture is wrong, but 
only that dropping the AQ term to deduce formula (JSJ) is not a valid assumption. 



3.2. Kerr black holes 

Quasinormal frequencies of Kerr black holes were first studied by Leaver using continued fraction techniques 4] 
and subsequently investigated by other authors |34j . 

A systematic exploration of the behaviour of the first few overtones was carried out only some years ago, using 
Leaver's continued fraction method, by Onozawa |l9). who found some rather odd features as the mode damping 
increases. For example, Detweiler |35| showed that the first few modes having I — m have vanishing imaginary part 
and real part equal to m as the black hole becomes extremal (in our units, as a — > 1/2). He also showed analytically 
that there can be infinite solutions to the Teukolsky equation having a = 1/2 and to = m, which led to the suspicion 
that all modes having / = m should "cluster" on the real axis at u> = m as the black hole becomes extremal. An 
interesting outcome of Onozawa's investigation was that, for given values of I and m = I, there is at least one mode 
frequency which does not tend to m in the extremal limit. Onozawa also found that, when n ~ 10, modes having 
negative m (which for the first few overtones show a tendency to decrease in frequency as the hole is spun up) show 
instead a tendency to "turn around" and increase their frequency as a increases, sometimes showing strange "loops" 
in the complex-w plane (see figures 3 and 4 in 0]). 

We have confirmed Onozawa's results, and extended them to moderately high n. However, since one has to solve 
simultaneously the angular and the radial continued fraction, the numerical problem turned out to be much more 
tricky than the search of highly-damped modes in the RN case. Comparing Nollert's technique and a standard 
summation of the continued fraction using Gautschi's algorithm H, we found out that it is much harder to achieve 
a stable numerical computation of modes for n ~ 50 or higher |36| . Because of these convergence problems, even 
using Nollert's method, we did not manage to push the numerical calculation to very large values of n. However, even 
moderately high values of n shed some light on what should be the asymptotic behaviour of Kerr QNM's. 

Let us first consider Kerr perturbations having m — 0. We have been able to compute quite a few moderately 
highly damped modes for I — 2 and I = 3, and in figure 01 we show two of these modes. These plots should be 
compared to the RN modes we have shown in figure ^ there is a similar "looping" behaviour, with the number of 
loops increasing as the damping of the mode increases. We notice that a similar looping behaviour has recently been 
found by Glampedakis and Andersson for scalar perturbations of Kerr black holes, using a different method 37]. 
Figure [S] and show the real and imaginary parts (respectively) of some Kerr modes with m = 0. It is useful to 
compare these plots with figures and [21 the behaviour is extremely similar, and asymptotically it can probably be 
described by some formula reminiscent of J()J|. Of course, such a behaviour is not even close to that predicted by 
formula JSJ. 

Hod recently used Onozawa's data, which to our knowledge are the only available published data for highly damped 
Kerr modes, to show that the results predicted by formula agree with the numerically computations in [19| within 
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~ 5%, at least when I = m. We repeated Onozawa's calculations, finding excellent agreement with his results, and 
then extended it to higher-order modes. We found that, as n increases, the formula conjectured by Hod does not 
seem to provide a good fit to the asymptotic modes. As shown in figure 03 the proposed formula disagrees quite 
badly in the small rotation rate regime even for the low-lying modes which were used for comparison by Hod. In any 
event, this is to be expected and does not contradict Hod's claim, since it is known from previous calculations in the 
Schwarzschild limit 4] that modes with n < 10 have real parts of the QNM frequencies which are not close enough to 
their asymptotic value. However, if formula (JjJJ really holds in the asymptotic limit, one would expect the agreement 
with numerical calculations to get better as the mode's imaginary part increases; on the contrary, figure^Jshows that, 
as n increases, the agreement gets worse, even for large rotation rates (which used to show rather good agreement for 
n ~ 10). So Hod's formula does not seem to provide an accurate fit to the asymptotic frequencies. We actually found 
that, as the mode order increase, modes having I = m = 2 are fitted extremely well (relative errors being of order 
~ 0.1% when n ~ 50) by the relation 

w££La = 2n + i^Tjin. (21) 

It is interesting to note that, although the real part tends to a different asymptotic value, the spacing in the imaginary 
parts does agree with the value conjectured by Hod in |l8j |. Indeed, Hod put forward his conjecture observing that in 
the Schwarzschild case the asymptotic QNM spacing is given by 2ttTh ■ The convergence to the indicated asymptotic 
behaviour is faster for large values of a, and generally the agreement between formula i|21|) and our numerics is good 
for values of a > 0.1. The index n appearing in formula (|21|) depends, of course, on the labelling convention we use 
to "count" QNM's. 

Notice also that the modes' imaginary part does not show the typical "shift" (—1/4) which is present in formula 
© ■ In our opinion, this is another hint that the Schwarzschild limit of highly damped Kerr modes should be taken 
with special care, and that order of limits issues may be relevant in the asymptotic regime. 

Preliminary calculations show that modes having I = 2 and to = 1, to = — 1 and m = — 2 show a more complicated 
behaviour. We are currently trying to improve our understanding of highly damped Kerr QNM's using both analytical 
and numerical techniques. Some of our results will be shown in a separate paper |38j |. 



4. CONCLUSIONS 



In this paper we have numerically investigated the asymptotic behaviour of QNM's for charged (RN) and ro- 
tating (Kerr) black holes. We have first confirmed Nollert's results and extended them to scalar perturbations of 
Schwarzschild black holes. Our numerics are consistent with 

In 3 k s [(s 2 -l)-3l(l + l)] 

UJR = ^ + 7^ ' ( } 

where, within our accuracy, fco = —0.0970 = 3fc_2- Recently the constant fc_2 has been determined analytically - see 
formula l|19l) - and our result may be useful to determine fco as well. 

More importantly, our results for charged and rotating black holes do not agree with the simple behaviour predicted 
by Hod's conjecture for the real part of the frequency, as given in formula JSJ. We have shown that both the real and 
the imaginary part of RN QNM's as functions of charge display an oscillating behaviour. The oscillations start at 
smaller values of the charge and get faster as n increases. We have compared our numerical results to the predictions 
of the asymptotic formula © derived in 0] and found that, in general, they agree extremely well, especially for 
large values of the mode index n. The formula derived by Motl and Neitzke only fails to reproduce our numerical 
results in the limit of small charge. Quite likely, their expression has a wrong limit in the Schwarzschild case (th e 
real part tending to ln5/47r instead of ln3/47r) because finite- ro corrections become relevant as Q tends to zero [l7| . 
A computation of higher-order corrections to the asymptotic RN formula will probably give a final answer on the 
reasons for the small-Q discrepancy we observe. If we trust the predictions of formula © , as the agreement between 
analytical and numerical results suggests to do, the asymptotic frequency for extremal RN black holes is the same as 
for Schwarzschild black holes. This implies that dropping the AQ term in the first law of black hole mechanics, as 
done by Hod to conjecture the validity of formula (jSJ, is presumably not a justified assumption. 

Highly damped gravitational quasinormal frequencies of rotating (Kerr) black holes proved more difficult to compute, 
even using Nollert's method, and show a more complicated behaviour. For m = 0, we observe spirals in the complex 
plane reminiscent of RN modes, confirming the behaviour found in (3^ for highly damped scalar perturbations. For 
/ = m = 2, the behaviour is completely different. The asymptotic behaviour of Kerr modes having I = m = 2 is very 
well fitted by formula 1)21(1 . that we rewrite here: 



.Kerr 



2VL + i2irT H n, 
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where f2 is the angular velocity of the black hole horizon and Th its temperature. The convergence to the limiting 
value is faster when a is large, and the formula can be see as an extension of formula @ to Kerr modes having 
I = m = 2. We think that a calculation of finite- n corrections to the asymptotics may help explain both the faster 
convergence rate at large a, and the apparent disa gree ment with formula @ in the limit a — > 0. A more extensive 
investigation of asymptotic Kerr QNM's is ongoing |38| . 
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FIG. 1: The top two panels show the behaviour of the n = 5 and n — 10 QNM frequencies in the complex ui plane. The 
n = 10 mode "spirals in" towards its value in the extremal charge limit; the number of spirals described by each mode increases 
roughly as the mode order n. The panels in the second row show how the n = 10 spiral "unwinds" as the angular index I is 
increased (in other words, the asymptotic behaviour sets in later for larger Z's). Finally, the bottom panel shows a high-Z mode 
trajectory "pointing" to its limit as the charge becomes extremal. In all cases, we have marked by an arrow the frequency 
corresponding to the Schwarzschild limit (Q = 0). 
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FIG. 2: Real part of the RN QNM frequencies as a function of charge for n = 5, 10, 30, 60, 5000, 10000, 100000. As the 
mode order increases the computation becomes more and more time consuming, the oscillations become faster, and a good 
numerical sampling is rather difficult to achieve; therefore in the last plot we use different symbols (small squares, circles and 
triangles) to display the actually computed points. For n = 5000, 10000, 100000 we also compare to the prediction of the 
analytic formula @ derived by Motl and Neitzke [Tfil ]. The oscillatory behaviour is reproduced extremely well by their formula, 
but the disagreement increases for small charge: formula © does not yield the correct Schwarzschild limit. 
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FIG. 3: Imaginary part of the RN QNM frequencies as a function of charge for n = 10, 30, 60, 5000. For n = 5000 we also 
display the actually computed points, and compare to the prediction of the analytic formula ©. As for the real part, the 
oscillations are reproduced extremely well, but the disagreement with our numerical data increases for small charge. 




FIG. 4: Trajectories of Kerr modes having m — in the complex-o; plane. The left panel corresponds to Z = 3, n — 15 and 
the right panel to Z = 3, n = 20. The number of spirals increases with the mode order, as in the RN case. We have marked by 
an arrow the point in the plane corresponding to the Schwarzschild limit. 
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FIG. 5: Real part of Kerr modes having m — as a function of a. Labels indicate the corresponding values of I and of the 
mode order n. 
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FIG. 6: Imaginary part of Kerr modes having m = as a function of a. The left panel corresponds to I — 3, n = 15 and the 
right panel to I = 3, n = 20. 




FIG. 7: Relative error in the formula for the asymptotic frequency conjectured by Hod. The plots shows E = [uIr — ujr)/ujr, 
where uir is defined in formula 10, for increasing values of the mode index n, namely n — 12, 30, 40, 50 and I = m = 2. As n 
grows, the relative error in the conjectured asymptotic formula tends to be larger and larger. 



